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We report the first observation of the Dalitz decay rj —> -ye + e~ , based on a data sample of 
1.31 billion J/ijj events collected with the BESIII detector. The rf mesons are produced via the 
J/tf —> 77/ decay process. The ratio T(?/ — > 7 e + e - )/T(?/ —► 77 ) is measured to be (2.13 ± 
0.09(stat.) ± 0.07(sys.)) x 1CP 2 . This corresponds to a branching fraction B(r/' —>■ 7 e + e _ ) = (4.69 ± 
0.20(stat.) ± 0.23(sys.)) x 10~ 4 . The transition form factor is extracted and different expressions 
are compared to the measured dependence on the e + e~ invariant mass. The results are consistent 
with the prediction of the Vector Meson Dominance model. 


PACS numbers: 13.40.Gp, 14.40.Be, 13.20.Gd, 13.40.Hq 
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I. INTRODUCTION 


Electromagnetic (EM) Dalitz decays of light pseu¬ 
doscalar mesons, P —> 7 l + l~ (P = 7r°, t], rf\ l = e,p), 
play an important role in revealing the structure of 
hadrons and the interaction mechanism between photons 
and hadrons [1]. If one assumes point-like particles, the 
decay rates can be exactly calculated by Quantum Elec¬ 
trodynamics (QED) [2]. Modifications to the QED decay 
rate due to the inner structure of the mesons are encoded 
in the transition form factor (TFF) F(q 2 ), where q is the 
momentum transferred to the lepton pair, and q 2 is the 
square of the invariant mass of the lepton pair. A recent 
summary and discussion of this subject can be found in 
Ref. [3]. 

The knowledge of the TFF is also important in stud¬ 
ies of the muon anomalous magnetic moment, a ^ = 
( g M — 2)/2, which is the most precise low-energy test 
of the Standard Model (SM) and an important probe 
for new physics. The theoretical uncertainty on the SM 
calculation of ap is dominated by hadronic corrections 
and therefore limited by the accuracy of their determina¬ 
tion [4]. In particular, the hadronic light-by-light (HLbL) 
scattering contribution to a „ includes two meson-photon- 
photon vertices that can be related to the form factors in 
P —> 77* —> 'ye + e~ decays [4], Thus, models describing 
these transitions should be tested as precisely as possi¬ 
ble to reduce the uncertainty in the SM prediction for 

(5m - 2)/2- 

In this work, the Dalitz decay rf —> 'ye+e is measured 
for the first time. The differential decay width, normal¬ 
ized to the radiative decay width 7' —> 77, is [1] 


dT(rf -a 7 l+l~) 
dq 2 T{rj' —> 77) 



1 - 


4m? 


1 + 


= [QED(g 2 )] x \F(q 2 )\ 2 , 




( 1 ) 


where m v / and mi are the masses of the rf meson and 
the lepton, respectively; a is the fine structure constant; 
and [QED(g 2 )] represents the calculable QED part for a 
point-like meson. The TFF, F{q 2 ) : which is described by 
phenomenological models, can be experimentally deter¬ 
mined from differences between the measured di-lepton 
invariant mass spectrum and the QED calculation. In the 
Vector Meson Dominance (VMD) model [5], it is assumed 
that interactions between virtual photon and hadrons are 
dominated by a superposition of neutral vector meson 
states. One commonly used expression for the multi-pole 
form factor is [6]: 


v 2 g Vl my — q z — iT v m v 


( 2 ) 


where N is a normalization constant ensuring that 
F( 0) = 1; V = p,u>,(j)\ my, Ty are the masses and 


widths of these vector mesons; and cy 7 y and gv-y are 
the corresponding coupling constants. 

The parameter to be experimentally determined is the 
slope of the form factor 6, which is related to the effective 
virtual vector meson mass A by 


b 


dF 

dq 2 



( 3 ) 


In experiments, the single-pole form factor is generally 
used to extract the slope of the form factor. For the case 
of the rf , the pole is expected to lie within the kinematic 
boundaries of the decay. The square of the form factor 
is described by 


|F(g 2 )| 2 


A 2 (A 2 +7 2 ) 
(A 2 — q 2 ) 2 + A 2 7 2 


( 4 ) 


where the parameters A and 7 correspond to the mass 
and width of the Breit-Wigner shape for the effective 
contributing vector meson. To a first approximation, A ss 
M p « 0.7 GeV and 7 « T p « 0.12 GeV. 

For the rf Dalitz decay, only the process rf -A 
has been observed and the slope of the form factor was 
measured to be £y = (1.7 ± 0.4) GeV -2 [1, 7]. To date, 
the process rf -A je + e~ has not been observed yet. The 
most stringent upper limit on the ratio of decay widths 
T(r]' 7e + e _ )/T(r/ —> 77) is 4.1 x 10 -2 at the 90% 

confidence level (CL) from the CLEO Collaboration [8], 
which is above the predicted value of (2.06±0.02) x 10 -2 
from the modified VMD model [9]. 

In the VMD model, the TFF slope is expected to 
be b v > = 1.45 GeV -2 [10, 11], while for chiral pertur¬ 
bation theory it is ty = 1.60 GeV -2 [12]. A recent 
calculation based on a dispersion integral gives iy = 
1.53^1 GeV" 2 [13]. 

We report the first observation of the rf -A ye+e 
decay and the extraction of the TFF. The source of the 
rf mesons are radiative J/rp —> 77/ decays in a sample 
of 1.31 billion J/ip events (2.25 x 10 8 events were taken 
in 2009 [14] and 1.09 x 10 9 in 2012) [15] collected by 
the BESIII [16] at the BEPCII e + e“ collider. The rf —> 
77 decay events in the same data sample are used for 
normalization. 


II. THE BESIII EXPERIMENT AND MONTE 
CARLO SIMULATION 

BEPCII is a double-ring multi-bunch e + e“ collider 
running in the tau-charm energy region. The BESIII 
detector, described in detail in Ref. [16], has a geomet¬ 
rical acceptance of 93% of 47 t solid angle. It consists of 
a drift chamber (MDC), a time-of-flight (TOF) system, 
and an electromagnetic calorimeter (EMC), all enclosed 
in a super-conducting solenoid with 1.0 T (0.9 T in 2012) 
magnetic field. The small-cell helium based MDC pro¬ 
vides the tracking of charged particle and ionization en¬ 
ergy loss ( dE/dx ) measurement. The single cell position 
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resolution is 130 jim and the transverse momentum reso¬ 
lution is 0.5% at 1 GeV/c. The TOF system for particle 
identification (PID) is made of plastic scintillators. It 
has 80 ps time resolution in the barrel, and 110 ps in 
the end caps. The EMC is made of 6240 Csl (Tl) crys¬ 
tals. The energy resolution is 2.5% in the barrel and 
5% in the end caps for 1.0 GeV photons. Outside the 
solenoid, a muon chamber system made of 1272 m 2 resis¬ 
tive plate chambers detects muon tracks with momenta 
greater than 0.5 GeV/c. 

The GEANT4-based [17] simulation software BOOST 
includes the description of geometry and material of the 
BESIII detector, the detector response and digitization 
models, and also tracks the detector running conditions 
and performance. A Monte Carlo (MC) simulated sample 
of 1.2 billion J/ip inclusive decays is used to study poten¬ 
tial backgrounds. The production of the J/■ip resonance 
is simulated by the MC event generator KKMC [18]; the 
known decay modes are generated by evtgen [19, 20] 
with branching fractions set at the world average val¬ 
ues [21], while unknown decays are generated by LUND- 
CHARM [22]. The evtgen package is used to generate 
J/ip —> 777 ', r]' —> 7 e + e~ and 7 ' —> 77 events. The decay 
J/ip —> 77 ]' is generated with an angular distribution of 
1 + cos 2 9 1 , where 9 1 is the radiative photon angle rela¬ 
tive to the positron beam direction in the J/ip rest frame. 
In generating r)' —> 7 /e + e~ , the TFF is parameterized by 
the multi-pole VMD model in Eq.(2) with the parameters 
taken from Ref. [1]. 

III. SIGNAL SELECTION: J/ip ->• 777', 77' ->• ye+e" 

Charged tracks are reconstructed from hits registered 
in the MDC. Only tracks with |cos0| < 0.93 are re¬ 
tained, where 9 is the polar angle with respect to the 
beam axis. The tracks are required to pass within 
10 cm of the center of the interaction region in the 
beam direction (Z axis) and within 1 cm in the plane 
perpendicular to the beam. Event candidates are re¬ 
quired to have two well reconstructed charged tracks with 
net charge zero. For electron identification, information 
from dE/dx and TOF is combined to compute probabil¬ 
ities for the electron (Prob(e)) and pion (Prob( 7 r)) hy¬ 
pothesis. To separate electrons from pions, we require 
Prob(e)/(Prob(e) + Prob( 7 r)) > 0.95. Final states with 
kaons cannot contribute to the background because of 
the limited phase space. 

Electromagnetic showers are reconstructed from clus¬ 
ters of energy deposits in the EMC. The photon candi¬ 
date showers must have a minimum energy of 25 MeV 
in the barrel region (| cos0| < 0.80) and 50 MeV in the 
end cap region (0.86 < |cos0| < 0.92). Showers in the 
region between the barrel and the end caps are poorly 
measured and excluded from the analysis. To exclude 
charged-particle induced activities, the showers are re¬ 
quired to be separated from the extrapolated positions 
of any charged track by at least 10°. In addition, cluster 


timing requirements are used to suppress electronic noise 
and unrelated energy deposits. 

A vertex fit is performed on the electron and positron 
tracks, and a loose y 2 requirement is applied to ensure 
that they come from a common vertex. To improve res¬ 
olution and reduce background, a four-constraint (4C) 
kinematic fit is performed to the 77 e + e _ hypothesis that 
constrains the total four-momentum of the detected par¬ 
ticles to be equal to the initial four-momentum of the 
colliding beams. For events with more than two pho¬ 
ton candidates, the combination with the smallest y| c ' s 
selected. Only events with y 2 c <100 are retained. 

For the J/ip xv'iV' —•* 7 e + e _ signal channel, 

the largest background comes from QED processes and 
J/ip —► e + e _ 77 decays. For these channels, the com¬ 
bination of the e + e _ with any final-state photon pro¬ 
duces a smooth M('ye + e~) distribution. The QED back¬ 
ground mainly comes from e + e _ —> e + e _ 77 and e + e _ —> 
37 events in which one 7 converts into an e + e~ pair. 
These are studied using a e + e _ collision data sample of 
2.92 fb -1 taken at y/s = 3.773 GeV [23], which is dom¬ 
inated by QED processes. For those processes, most of 
the photons have low energy and are at small angles rela¬ 
tive to the incoming electron or positron beam directions. 
To reduce this background, the energy of the low-energy 
photon is further required to be higher than 200 MeV, 
and the angle between the photon and the electron or 
positron initial direction in the final states is required to 
be larger than 10 °. 

The primary peaking background comes from the de¬ 
cay J/ip —> 777 ', 77 ' —> 77 followed by a 7 conversion in the 
material in front of the MDC, including the beam pipe 
and the inner wall of the MDC. The distance from the re¬ 
constructed vertex point of the electron-positron pair to 

the z axis, defined as 5 xy = R x + R y , is used to distin¬ 
guish 7 conversion events from signal events [24], where 
R x and R y are the distances in the x and y directions, 
respectively. A scatter plot of R y versus R x is shown 
in Fig. 1 (a) for MC-simulated J/ip —»■ 777', 77' —> 77 
decays, in which one of the photons undergoes conver¬ 
sion to an e + e _ pair. As indicated in Fig. 1 (a), the 
inner circle matches the position of the beam pipe, while 
the outer circle corresponds to the position of the inner 
wall of the MDC. Figure 1 (b) shows the 6 xy distribu¬ 
tions for the MC simulated J/ip —> 777' ,77' —> ye+e - , 
J/ip 777', 77' —> 77 events, together with the selected 
data events and events from the v/ mass sideband. The 
two peaks above 2.0 cm correspond to the photon con¬ 
version of the 7 from J/ip —> 777' ,77' —> 77 events, while 
the events near 5 xy = 0 cm originate from the inter¬ 
action point. We require S xy < 2 cm to suppress the 
photon-conversion background, which retains about 80% 
of the signal events while removing about 98% of the 
photon-conversion events. After all selections, the nor¬ 
malized number of expected peaking background events 
from J/ip —>■ 777', 77' — > 77 is 42.7 ± 8.0, where the error 
is dominantly from the difference in selection efficiencies 
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for the 7 conversion events between data and MC. 

Another possible source of peaking background is 
J/ip —> '■ 77 /,?/ —> 77 r + 7 r _ , where the two pions are 
misidentified as an e + e~ pair. The momenta of these 
pions are most below 200 MeV/c. An exclusive MC sam¬ 
ple that includes coherent contributions from p, w, and 
the box anomaly in the decay rj' —> 'yTT + n~ [25] is used 
to study this background. We find that the kinematic fit 
to the electron-positron hypothesis shifts the spectrum 
away from the rj' mass and, thus, the resulting M( 7 e + e _ ) 
distribution does not peak at the r\ mass value. The nor¬ 
malized number of events from this background source af¬ 
ter all selections is 9.7±0.4, which is negligible compared 
to the non-peaking background from e + e _ -A e + e~"f"f. 



-10 -5 0 5 10 0 2 4 6 8 10 


Rx (cm) 8 xy (cm) 

Figure 1 . Electron-positron vertex position distribution: (a) 
scatter plot of R y versus R x for MC-simulated J/ip —> 7?/, 
7/ -A 77 events, (b) 8 xy distributions. The (black) crosses 
are data. The (red) dashed line shows the MC-simulated 
J/ip — y 777', r\ — > 7e + e“ signal events. The (orange) dotted- 
dashed histogram shows the background from 7-conversion 
events. The (green) shaded area is estimated from the r/' mass 
sideband. The (blue) line is the sum of MC and the sideband 
estimate. In (b), the solid arrow indicates the requirement on 

8xy • 


The combination of 7 e + e _ with invariant mass clos¬ 
est to m v ' is taken to reconstruct the 77 '. The result¬ 
ing Mpye + e~) distribution after the selection criteria 
is shown in Fig. 2 and exhibits a clear peak at the 7 / 
mass. An unbinned extended maximum likelihood (ML) 
fit is performed to determine the signal yield. The signal 
probability density function (PDF) is represented by the 
signal MC shape. The shape for the non-peaking back¬ 
ground is described by a first-order Chebychev polyno¬ 
mial. The background yield and its PDF parameters are 
allowed to vary in the fit. The peaking background from 
the 7 -conversion of J/ip — > 777 ', —> 77 decay is obtained 
from the MC-simulated shape with the yield fixed as de¬ 
scribed before. The fitting range is 0.85—1.05 GeV/c 2 . 
The net signal yield and the detection efficiency are sum¬ 
marized in Table I. 


Table I. Number of observed signal events, N v ,^ e + e - 
(A^/_> 77 ), and detection efficiency, e v r^ e + e - (e^'_> 77 ) for 
J/ip -A 777', r / 1 -A -ye + e~ (J/ip -A 777', 77' “A 77 )- The un¬ 
certainties are statistical only. 



77' -a 7e+e 

if -A 77 

^r]'— l7e+e“ (-^77'—^77) 

864 ± 36 

70846 ± 292 

^77'—l7e+e“ (^77'—^77) 

24.5% 

42.8% 



lected signal events. The (black) crosses are the data, the 
(red) dashed line represents the signal, the (green) dot-dashed 
curve shows the non-peaking background shapes, the (orange) 
shaded component is the shape of the J/ip -A 777', 77' —> 77 
peaking background events. The total fit result is shown as 
the (blue) solid line. 


IV. NORMALIZATION CHANNEL: 

J/ip — > 777', 77' — y 77 

The decay J/ip -A 777', 77' -a 77 is studied using the 
same data set, and we quote our result in terms of the 
ratio r(7?' -A 7e + e _ )/r(77' -a 77). In this ratio the un¬ 
certainties due to the total number of J/ip events and 
the branching fraction for J/ip -A 777' cancel, and the 
uncertainty due to the photon detection efficiency par¬ 
tially cancels. 

Events with zero charged particles and at least three 
photon candidates are selected with the same require¬ 
ments that are used for the signal events. A 4C kine¬ 
matic fit is performed to the J/ip -a 777 hypothesis. For 
events with more than three photons, the combination 
with the smallest x\c selected. The xlc required 
to be less than 100. The two photon combination with 
invariant mass M (77) closest to m v / is taken as from the 
77 ' decay. 

Detailed MC studies indicate that no peaking back¬ 
ground remains after all the selection criteria. The non¬ 
peaking background mainly comes from the continuum 
process e + e~ -A 777 and J/ip -A 77r°7r° decays. The 
latter source involves intermediate states such as the 
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/o(1500), /o(1710), / o (2020), / 2 (1270), / 4 (2050). Because 
the T]' decays isotropically, the angular distribution of 
photons from the 77 ' decays is flat in cos 0 decay ; where 
^decay is the angle of the decay photon in the 77 ' helic- 
ity frame. In contrast, background events from QED 
continuum processes and J/ij; —> "fn 0 n° decays tend to 
accumulate near cos 9decay = ± 1. We suppress these non¬ 
peaking backgrounds by requiring | cos 0 <jecay| < 0 . 8 . 

The M ( 77 ) distribution for events that survive the se¬ 
lection requirements is shown in Fig. 3. An unbinned ML 
fit is performed to obtain the yield of J/i8 —> 777 ', t?' —> 
77 . The PDF used to represent the signal is taken from 
the MC, and the PDF for the non-peaking background 
is a first-order Chebychev polynomial with coefficients 
determined from the fit. The resulting signal yield and 
the MC-determined detection efficiency are summarized 
in Table I. 



Figure 3 . Invariant mass distribution of M( 77) for the se¬ 
lected events in the normalization channel. The (black) 
crosses are data, the (red) dashed curves represents the 
v' “A 77 signal, the (green) dot-dashed curve shows the non¬ 
peaking background. The fit result is shown as the (blue) 
solid curve. 


V. SYSTEMATIC UNCERTAINTIES IN THE 
RELATIVE DECAY WIDTH 

Table II lists all sources of systematic uncertainties 
associated with the measurement of the ratio T(r)' —> 
7 e + e _ )/r(? 7 ' —> 77 ). Most systematic uncertainties are 
determined from comparisons of low-background, high- 
statistics data samples with results from MC simulations. 

The electron and positron tracking and PID efficien¬ 
cies are determined using a sample of radiative Bhabha 
e + e _ -A "fe + e~ (including J/if) -a 7 e + e _ ) events col¬ 
lected at the J/t/> energy. Differences in tracking and 
PID efficiencies between data and MC simulation are de¬ 
termined for every bin of a two-dimensional distribution 
of the momentum versus polar angle of the lepton tracks. 
These are used to determine an overall weighted differ¬ 
ence per track of (1.1 ± 0.3)% for the tracking efficiency 
and (1.9±0.3)% for the PID efficiency. The MC efficiency 
is corrected for these differences, and the uncertainties 


Table II. Summary of relative systematic uncertainties for the 

determination of the ratio r C ~) 7e e ) _ The i as t row j s the 

"(A—> 77 ) 

uncertainty associated with the J/tj) — > 777', v[ -a 77 normal¬ 
ization sample. 


Sources 

Systematic uncertainties(in %) 

MDC tracking 

0.6 

PID 

0.6 

Photon detection 

1.0 

Veto of gamma conversion 

1.8 

4 C kinematic fit 

1.0 

Form factor uncertainty 

1.6 

Fit range & Bkg shape 

0.9 

Uncertainty of A^/_> 77 

1.2 

Total 

3.3 


of the correction coefficients are assigned as the system¬ 
atic uncertainties associated with the lepton tracking and 
PID efficiencies. 

The systematic uncertainty associated with the photon 
detection efficiency is studied using three different meth¬ 
ods, as described in Ref. [26]. The three methods provide 
consistent results for the photon efficiency uncertainty to 
be 1% per photon. Because the systematic uncertainty 
from the radiative photon and one photon from the rj' 
cancel in the ratio, the total systematic uncertainty from 
photon detection is 1 %. 

In the analysis, the peaking background from J/iji -A 
777', 77 ' -A 77 7 -conversion events is suppressed by the 
requirement 8 xy < 2 cm. To estimate the systematic un¬ 
certainty associated with this requirement, we use a sam¬ 
ple of J/t/j -a 7r + 7r _ 7r°, 7 t° -a 7 e + e _ that includes both 
7T° Dalitz decays and 7 r° —> 77 decays with one photon 
externally converted to an electron-positron pair. The 
data-MC difference of 1.8% for these events is considered 
as the systematic uncertainty for our 7 -conversion veto 
requirement on 5 xy . 

A systematic uncertainty associated with the kine¬ 
matic fit will occur if the track-helix parameters for data 
and MC simulated events are not consistent. Follow¬ 
ing the procedure described in Ref. [27], we use the 
J/ip -A 7r + 7r _ 7r°, 7T° -A 7e + e _ decay as a control sample 
to extract the correction factors from the pull distribu¬ 
tions of the track helix parameters. The 1% difference 
between the efficiencies with and without helix parame¬ 
ter corrections is taken as the systematic uncertainty. 

To estimate the systematic uncertainty due to the effi¬ 
ciency dependence on the choice of form-factor param- 
eterizations, signal MC events are also generated us¬ 
ing a single-pole VMD model, shown in Eq. (4), with 
A = (0.79 ±0.05) GeV and 7 = (0.13 ±0.06) GeV, which 
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are taken from the fitted results described below in Sec¬ 
tion VII. The relative difference in the detection efficiency 
compared to that of the multi-pole model is taken as the 
uncertainty associated with the form-factor parameteri¬ 
zation. 

In the fit to the 7 e + e“ mass distribution, the signal 
PDF is fixed to the signal MC shape. An alternative 
fit is performed by using a convolution of a MC signal 
shape with a Gaussian function that is used to describe 
the MC-data difference due to the resolution. The fitted 
width of the Gaussian is (0.39 ± 0.19) MeV, and the fit 
yields 863.8 ± 36.0 signal events. The difference from 
the nominal fit is negligible. Finally, the uncertainty due 
to the non-peaking background shape is estimated by 
varying the PDF shape and fitting range in the ML fit. 
The changes in yields for these variations give systematic 
uncertainties due to these backgrounds. 

The systematic uncertainty in the measurement of 
J/ij) 77 /, rj' — > 77 associated with the uncertainty from 
the kinematic fit is estimated using a control sample of 
e + e - —► 777 at 3.650 GeV [28] and found to be less 
than 1%. The uncertainty for this channel due to back¬ 
ground is estimated to be less than 0.3% from variations 
in the PDF shape and fitting range. The uncertainty 
from the requirement |cos0<j e cay| < 0.8 is 0.4%. When 
combined with the 0.4% statistical uncertainty, the total 
uncertainty associated with A ^/_> 77 is 1 . 2 %. 

Assuming all systematic uncertainties in Table II are 
independent, the total systematic uncertainty, obtained 
from their quadratic sum, is 3.3%. 


VI. RELATIVE DECAY WIDTH 


The ratio T(r]' -A je + e )/T (7/ — > 77) is determined 
using the following formula: 

T(7 17 e + e ) _ V^/_^ 7e + e — e^'_^ 77 , 

r(v' -> 77 ) N v *11 yye+e~ 

where A r I? /_ ) . 7e + e - (AT^/_^ 77 ) and e ?; /_ ) , 7e + e - (e T7 /_ s . 77 ) are 
the number of observed signal events and the detec¬ 
tion efficiency, respectively, for J/ip —> ^jrf ,r}' -a 7e + e - 
(J/ip —► 77 /, 7/ -a 77) decays, as listed in Table I. The 
result is 

r (7? ->7e+e ) = ( 2 . 13 ± 0 . 09 (stat.)± 0 . 07 (sys.)) x 10" 2 . 
IW -A 77) 

( 6) 

Using the rf —> 77 branching fraction value listed in 
PDG [21], we obtain the first measurement of the 77 ' —> 
7e + e - branching fraction of 

23(77' -a 7e+e") = ( 4 . 69 ± 0 . 20 (stat.)± 0 . 23 (sys.)) x 10 -4 . 

( 7 ) 


VII. FORM FACTOR MEASUREMENT 


The TFF is extracted from the bin-by-bin efficiency 
corrected signal yields for eight different M(e + e~) bins. 
The bin widths are all chosen to be 0.1 GeV/c 2 . Since 
this is much wider than the M{e + e~) resolution, which 
is 5~6 MeV/c 2 depending on M(e + e~), no unfolding is 
needed. The signal yield in each M(e + e~) bin i is ob¬ 
tained by performing bin-by-bin fits to the M('ye + e~) 
mass distributions using the fitting procedure described 
in Section III. The peaking background from the J/ip -A 
777', r\' — > 77 only exists in the first bin, and the yield is 
fixed to the normalized number in the bin. The fit re¬ 
sults are shown in Fig. 4. The fitted (n“ bs ) and efficiency- 
corrected signal yields (n) :orr ) for each M,(e + e - ) bin are 
summarized in Table III. Figure 5 shows the efficiency- 
corrected signal yields versus M{e + e~ ) with the QED 
shape superimposed for comparison. The discrepancy 
between QED and data, which reflects the TFF, is evi¬ 
dent in the high M{e + e~ ) region. 

The systematic uncertainties on n£ orr include the un¬ 
certainty from the MDC tracking efficiency, PID, photon 
detection, kinematic fit, veto of gamma conversion, back¬ 
ground description and signal shape; they are the same 
as those described in Section V. 

The partial ratio 77 = Ar(7/ -A "je + e~)i/T{ri' -A 77) 
for each given M(e + e~ ) bin i, is defined as 

= AT(t/ -a 7e + e~)j = n^ orr e 7j /^ 77 
i - T(rf -a 77) ify- 77 ’ 1 J 


where Ar(?/ -a 7 e + e - )i is the integrated rate in each 
M(e + e~) interval. 

The result for |F | 2 in each M[e + e ~) bin is obtained by 
dividing the value 77 by the integrated QED predication 
in each M{e + e~) interval (see Eq. (1)). The values of 
|E | 2 for each M(e + e~) bin are summarized in Table IV. 

A variety of models have been traditionally used to 
parameterize the TFF. The most common one, based on 
VMD [5], uses only the first term in the dispersion re¬ 
lation. In this single pole model, the TFF is given by 
Eq. (4). The results of a least-squares fit with the single 
pole model is shown in Fig. 6 ; the parameters of the form 
factors are determined to be A r( / = (0.79 ± 0.05) GeV, 
7 y = (0.13 ± 0.06) GeV. From the fitted value of the 
parameter A^/, the slope of the form factor is obtained 
to be (1.60 ± 0.19) GeV -2 , in agreement with the re¬ 
sult fy = (1.7 ± 0.4) GeV -2 obtained in the process of 
11 -A [1] - 

To test the robustness of the slope extracted from the 
simple pole model, we also fit the data below 0.5 GeV/c 2 
using the single pole Ansatz used in lighter meson studies: 


F(q 2 ) 


1 

(l-9 2 /A 2 )’ 


(9) 


The parameterization diverges at M(e + e ~) = A and, 
therefore, can not be used for the whole kinematic region. 
The result of this fit is shown in Fig. 7. The slope of the 
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Figure 4. Results from bin-by-bin fits to the M('ye + e~) distributions for different M(e + e~) bins. The (black) crosses are data, 
the (red) dashed curves represent the signal, the (green) dot-dashed curves show the non-peaking backgrounds, the (orange) 
shaded component for the M(e + e~) < 100 MeV/c 2 bin is the shape of the peaking background from J/ij) —> 7 rj , j \ ' —> 77. The 
total fit results are shown as (blue) solid curves. 


Table III. Fitted (n° bs ) and efficiency-corrected (nf orr ) signal yields for the eight M(e + e ) bins, and ratios ( 77 ). The uncer¬ 
tainties are statistical only. 


M(e + e~) (GeV/c 2 ) 

[0.0, 0.1] 

[0.1, 0.2] 

[0.2, 0.3] 

[0.3, 0.4] 

n° bs 

545 ± 27 

86.5 ±10.7 

62.1 ±9.8 

45.6 ±9.7 

nf rr 

2380 ± 120 

368 ± 46 

194 ± 31 

128 ± 27 

ri (KT 2 ) 

1.44 ±0.07 

0.22 ±0.03 

0.12 ±0.02 

0.08 ±0.02 

M(e+e~) (GeV/c 2 ) 

[0.4, 0.5] 

[0.5, 0.6] 

[0.6, 0.7] 

[0.7, 0.8] 

n° bs 

45.4 ±9.0 

29.9 ±8.0 

28.0 ±7.8 

25.0 ±6.9 

nf rr 

135 ± 27 

93.3 ±25.0 

96.2 ± 26.8 

109 ± 30 

n (10“ 2 ) 

0.08 ± 0.02 

0.06 ±0.02 

0.06 ± 0.02 

0.07 ±0.02 
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Table IV. Values of |F| 2 in each M(e + e ) bin, where the first uncertainties are statistical and the second ones systematic. 


M(e+e—) (GeV/e 2 ) 

[0.0, 0.1] 

[0.1, 0.2] 

[0.2, 0.3] 

[0.3, 0.4] 

|F| 2 

1.05 ±0.05 ±0.03 

1.12 ±0.14 ±0.04 

1.16 ±0.18 ±0.05 

1.33 ±0.28 ±0.05 

M(e+e~) (GeV/e 2 ) 

[0.4, 0.5] 

[0.5, 0.6] 

[0.6, 0.7] 

[0.7, 0.8] 

|F| 2 

2.48 ± 0.49 ± 0.25 

3.30 ±0.88 ±0.31 

7.66 ±2.13 ±0.89 

26.6 ±7.3 ±1.9 



M(e + e ) (GeV/c 2 ) 


Figure 5. Efficiency-corrected signal yields ri; orr versus 
M(e + e~). The (black) crosses are data and the (gray) shaded 
histogram indicates the point-like QED result. 

form factor is determined to be b v > = (1.58±0.34) GeV -2 , 
which is in good agreement with the result of (1.60 ± 
0.19) GeV " 2 using Eq. (4). 

The quadratic difference between the uncertainties of 
the parameters with only statistical errors used in the 
fits and the uncertainties of the parameters with com¬ 
bined statistical and systematic errors used in the fits 
is taken as the systematic uncertainty on the param¬ 
eters. The resulting parameters in Eq. (4) are deter¬ 
mined to be A v r = (0.79 ± 0.04(stat.) ± 0.02(sys.)) GeV, 
7 y = (0.13 ± 0.06(stat.) ± 0.03(sys.)) GeV, respectively. 


VIII. SUMMARY 

In summary, with a sample of 1.31 billion J/tp events 
collected in the BESIII detector, we have made the 
first measurement of the EM Dalitz decay process rf —> 
"fe + e~ and measure the ratio T(r/ —> 7 e + e - )/r(?/ —> 
77 ) = (2.13 ± 0.09(stat.) ± 0.07(sys.)) x 10 -2 . Us¬ 
ing the PDG value for the rf —> 77 branching frac¬ 
tion [21], we determine B{rf -a 7 e + e - ) = (4.69 ± 
0.20(stat.) ± 0.23(sys.)) x 10 -4 . We present measure¬ 
ments of the TFF as a function of M(e + e~). Our 
TFF results can be described with a single pole pa¬ 
rameterization Eq. (4), with mass and width parame¬ 


ters of Ay = (0.79 ± 0.04(stat.) ± 0.02(sys.)) GeV, and 
7 y = (0.13 ± 0.06(stat.) ±0.03(sys.)) GeV, respectively. 
The slope of the TFF corresponds to (1.60±0.17(stat.)± 
0.08(sys.)) GeV -2 and agrees within errors with the 
VMD model predictions. The uncertainty of the rf tran¬ 
sition form factor slope matches the best determination 
in the space-like region from the CELLO collaboration 
&y = (1.60 ± 0.16) GeV -2 [29], and improves the pre¬ 
vious determination of the slope in the time-like region 
= (1.7 ± 0.4) GeV -2 [1, 7]. The rf form factor is 
determined by both universal 7 r + 7 r - rescattering and a 
reaction specific part, with the latter contributing about 
20% to the form factor slope [13]. Therefore our result 
is sensitive specifically to the rf internal EM structure. 
In addition, the decay rf -A ye+e - is closely related to 
rf —> 77 r + 7 r - , and in particular the transition form factor 
could be predicted from the invariant mass distribution of 
the two pions and the branching ratio of the rf —> 77 r + 7 r - 
decay in a model independent way using a dispersive inte¬ 
gral. Also, the knowledge of the TFF is useful for studies 
of the HLbL scattering contribution to the muon anoma¬ 
lous magnetic moment, a M = ( g^ — 2)/2 [4]. 
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Figure 6. Fit to the single pole form factor IF 1 ! 2 using Eq. (4). 
The (black) crosses are data, where the statistical and system¬ 
atic uncertainties are combined, the (blue) solid curve shows 
the fit results. The (gray) dotted line shows to the point-like 
case ( i.e. with |F| 2 = 1) for comparison. 
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Figure 7. Determination of the form factor slope by fitting to 
|F| 2 using the single pole form factor of Eq. (9). The (black) 
crosses are data, where the uncertainties are the combined 
statistical and systematic uncertainties, the (blue) solid curve 
shows the fit result. The (gray) dotted line corresponds to the 
point-like particle case (with |F| 2 = 1). 
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